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EFFECT OF POLLUTANT GASES ON OZONE PRODUCTION 


BY SIMULATED SOLAR RADIATION 
by E. L. Wong and David A. Bittker 
Lewis Research Center 

INTRODUCTION 

The present investigation provides some information related to the 
problem of ozone destruction by nitric oxide, NO, emission from high 
flying vehicles. According to Johnston (1971) and Westenberg (1972), 

NO emission could catalytically destroy ozone, Og, in the upper atmo- 
sphere to such an extent that additional ultraviolet, u. v. , radiation 
would reach the earth with possible harmful effects. 

F. E. Belles of this Center suggested that if Johnston and Westen- 
berg had included certain photochemical smog reactions in their NO-Og 
destruction scheme, NO destruction of Og might be reduced. Rapid NO 
conversion to nitrogen dioxide, NOg, followed by production of Og, has 
been observed in smog chamber experiments and Los Angeles smog 
(Leighton (1961)). This rapid NO conversion process to NOg followed 
by Og production is now believed to be a chain reaction involving un- 
saturated hydrocarbon and/or carbon monoxide, CO, oxidation (Leighton 
(1961), Heicklen et al. (1969), Westberg and Cohen (1971)). Thus, if 
engine exhausts have insufficient amounts of unsaturated hydrocarbon 
and CO to counteract the destruction of Og by the NO emitted, additional 
appropriate amounts of unsaturated hydrocarbon and/or CO could be 
added to engine exhaust gases. 

Based on Belles 7 recommendation, experiments were conducted 
using simulated solar radiation in a chamber with controlled atmosphere. 
Pressure was near 1 atmosphere in these tests. Tests were made to 
evaluate NO destruction of Og in the chamber and to determine if addition 
of CO and HgO could reduce NO destruction of Og. In conjunction with 
these experiments, the General Chemical Kinetics Program (Bittker and 
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Scullin (1972)) developed at the Lewis Research Center was modified 
and extended to permit a detailed kinetics data analysis. 



The recombination reaction O + O + M — Og + M also occurs, but is 
of negligible importance at atmospheric and stratospheric conditions. 
In this discussion we neglect the formation of excited state species in 
the photolysis reactions. If it is assumed that oxygen atom concentra- 
tion rapidly reaches a steady state value, then one can derive a differ- 
ential equation for Og production in which all terms containing are 
negligibly small. This equation can be integrated to give the following 
expression 




(i) 
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where 

t = reaction time 
[Og] = molar Og concentration 
[M] = total molar concentration 
v 2 = k 2 [0 2 ][M] 

In the expression e~ x , the exponent x is smaller than 0. 1 for reaction 
times up to about 1000 hours, so that the expansion 

e" x = 1 - x 

can be used. When this is done, the above equation becomes 


[O s ]= 2J 1 [0 2 ]t (2) 

Therefore, when pure air is irradiated, Og concentration should rise 
linearly for a long time. The rate of rise is determined only by the 
rate of photolysis of Og to form oxygen atoms. 

When dry air is irradiated in the presence of oxides of nitrogen, 
the following additional reactions occur which tend to destroy Og: 

k 5 

NO + Og-^NOg + 0 2 (V) 

J 6 

hv + N0 2 — NO + O (VI) 

k 7 

o + no 2 -^ no + o 2 (VII) 

k 8 

no 2 + O,— »• NO, + o 2 (VIII) 



4 


2NO + Og— *• 2N0 2 (IX) 

k 10 

NO + O + M — » NO, + M (X) 

J 11 

hi/ + NO, — NO + 0 2 (XI) 

Johnston (1971) and Westenberg (1972) present steady-state mathemati- 
cal analyses with the above mechanism. They show that the final 
concentration level is indeed much lower than it is when oxides of 
nitrogen are not present. However, in assessing the effect of pollutants 
injected into the atmosphere, it is necessary to obtain the time variation 
of Og concentration from the instant of injection. To obtain such concen- 
tration profiles the chemical mechanism must be numerically integrated. 
To perform this kind of computation, a general chemical kinetics program 
was developed several years ago at the Lewis Research Center (Bittker 
and Scullin (1972)). It has recently been expanded and improved to 
allow the inclusion of more types of reactions, including photochemistry 
and to make the integration method more efficient. This program_has 
been used extensively in the present work to analyze and interpret the 
experimental results. 

As mentioned previously the presence of H 2 0 and CO in jet engine 
exhaust gases gives rise to a series of reactions which may counteract 
the destruction of ozone (Heicklen et al. (1969)). 


k 12 


CO + OH — * C0 2 + H 

(XII) 

k 13 

H + 0 2 + M >H0 2 + M 

(xra) 


NO + H0 2 


k 


14 

-> N0 2 + OH 


(XIV) 
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H + NO, 


NO + OH 


(XV) 


Reactions XII through XIV are quite fast and provide a strong competi- 
tive path for the conversion of NO to NO 2 without the destruction of O 3 . 
Reaction XV, which could regenerate NO, was found to be unimportant 
under our experimental conditions. One must also consider the possi- 
ble destruction of O 3 in the presence of HgO alone, according to the 
chain mechanism: 


OH + O, 


ho 2 + 0 2 


(XVI) 


ho 2 + 0 3 


OH + 2 O, 


(XVII) 


The net effect of these last reactions is the conversion of two molecules 
of O 3 to three molecules of oxygen. In fact, H 2 0 and not NO was origi- 
nally thought to be the main destroyer of O 3 in engine exhaust gases 
(Kellogg (1970)). The most recent information on the values of k^g 
and k^ indicates that they are much smaller than originally estimated 
(Garvin and Hampson (1974) and DeMore (1973)). The destruction of O 3 
by NO is now considered much more important than its destruction by 
the water mechanism. Nevertheless the mechanism of reactions XVI 
and XVII cannot be ignored when studying atmospheric chemistry. The 
effect of water on O 3 production has been studied briefly in the present 
work. 

The formation of nitric and nitrous acids, HNO^ and HN0 2 , has 
been suggested by many people as a means of removing NO from the 
stratosphere and thus reducing its destructive effect on ozone (Crutzen 
(1972), DeMore (1973), Dutsch (1972)). The possible formation of these 
species in the gas phase during the present work was checked by theo- 
retical computations using these additional reactions: 
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k 18 

NO + OH + M — * HN0 2 + M (XVIII) 


k 19 

oh + hno 2 — »■ h 2 o + no 2 (XIX) 


k 20 

N0 2 + OH + M — * HNO, + M (XX) 


k 21 

OH + HNO, — *■ HgO + NO, (XXI) 


APPARATUS AND MATERIALS 

The main apparatus for this experimental work consists of a cham- 
ber. a solar simulator, and instrumentation for monitoring and NO x 
(NO and NO 2 ) concentrations. 

Reaction Chamber 

An existing chamber at this Center was modified to permit addition 
of pollutant gases and gas sampling (Fig. 1). The chamber has a volume 
of about 650 liters (23 cu ft) with approximate dimensions of 82 cm 
(32^ in. ) inside diameter, and length of 122 cm (48 in. ). Radiation 
enters the chamber through a 30.5 cm (12 in. ) diameter quartz window 
at one end. In front of this window is a water cooled shutter so that 
radiation can be turned TT on” or "off" instantly. This chamber can be 
evacuated to the low 10"° torr range using a 30. 5 cm (12 in. dia. ) oil 
diffusion pump fitted with a liquid nitrogen cold trap. The inside cham- 
ber surface is stainless steel. The radiation, after passage through the 
30.5 cm window, is a slightly diverging beam inside the chamber, as 
shown in Fig. 1. At the opposite end of the chamber the radiation de- 
scribes a light circle of about 40. 6 cm (16 in. ) in diameter. This indi- 
cates that about one-fifth of the effective chamber volume was exposed 
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to the radiation. Two Teflon sampling probes were positioned to sampl< 
gas in the center of the light beam in the chamber (Fig. 1). 

Solar Simulator 

Simulated solar radiation was provided by a commercial solar simi 
lator. With the apparatus shown in Fig. 1, the amount of Og produced 
in the chamber when dry clean air is irradiated, depends to a large ex- 
tent on the type of xenon arc lamp used. Three types of 4. 2 kW xenon 
high pressure arc lamps - I. T. T. , Ushio, and Osram - compatible 
with our solar simulator were tested. The peak Og production in pphm 
was as follows: 

pphm 

I.T.T. 20-30 

Ushio 40-60 

Osram 120-140 

Due to its much higher Og production, an Osram lamp was used through 
out this study. According to the manufacturer such lamps will provide 
u. v. radiation beginning at about 200 nm. In the visible range its radia 
tion is comparable to that emitted by the sun. Ultraviolet transmission 
through the quartz window was measured only as low as 250 nm, but it 
is possible that some radiation as low as 200 nm is also present in the 
transmitted light. 

Analytical Instruments 

For this work two types of commercially available chemilumines- 
cence analyzers were used to monitor Og and NO x (NO + NOg). Both 
instruments can only be used to sample at or near atmospheric pres- 
sure condition (low pressure limit ~700 torr). Because of this limita- 
tion the work reported here was obtained for the near atmospheric 
pressure condition only. 
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Test Gases 

Clean dry air used in the irradiation experiment was commercially 
available tanked ultrapur e clean ("UPC") air. This "UPC" air 
generally has the following typical stated analysis for contaminants in 
ppm by volume: 

COg 0.5 NgO (nitrous oxide) 0. 1 

CO 1.0 Total hydrocarbon 0. 1 

CH 4 (methane) 1.0 H 2 0 <1.0 

Desired amounts of were added to "UPC" air by saturating 
the air with distilled HgO at some pre-selected low temperature 
just prior to its entry into the chamber. 

Two other tanked gases were used. A mixture of NO (18 ppm) in 
nitrogen was used to calibrate the NO analyzer and to introduce NO gas 

A 

into the chamber. CO had a stated purity of 99.5 percent, and mass 
spectrometer analysis indicated that the 0. 5 percent impurity in CO was 
mainly C0 2 . 


EXPERIMENTAL RESULTS 
Evaluation of Apparatus for Ozone Production 

— 6 

For the O^ production tests the tank is first pumped to the low 10 
torr range for about a 20 hour period. Then "UPC" air is used to fill 
the tank to near atmospheric pressure. At time "zero, " irradiation of 
this air is begun by removal of the shutter over the quartz window. The 
Og concentration in the tank is monitored as a function of time. Two 
production tests were made in succession on the same day, and these 
are shown in Fig. 2 where concentration in pphm is plotted against 
irradiation time in minutes. During these tests the chamber surface had 
not been exposed to any pollutant gases such as CO, NO x , or H 2 0. For 
the two runs it is interesting to note that production for the second 
run was much more regular. Results of these and other runs suggest 
that a wall conditioning process was occurring. After this conditioning 
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process, Og production rate was less erratic. In the second run, peak 
steady state Og concentration was about 115 pphm in an irradiation time 
of about 40 minutes. Also on the same figure is shown the ozone decay 
(dark circled data) when irradiation is stopped. For the two runs the 
decay was an exponential one. Decay rates for the first and second run 
were about 0.011 min”* and 0.015 min~\ respectively. These results 
indicate that the Og decay may not be due entirely to wall reaction alone, 
otherwise one would expect the decay rate for the second run (after wall 
conditioning) would be less than that for the first run. 

Effect of Addition of NO to M UPC M Air on Og Production 

The effect of adding varying amounts of NO to "UPC" air on Og 
production is shown in Figs. 3, 4, and 5. For purpose of comparison, 
an Og production reference curve for air alone is shown in the 
three figures. This reference curve is an average of some Og 
production data for air alone obtained just before and after these tests 
involving addition of NO were made. Figure 3 shows that when 
air containing an initial concentration of 30 pphm NO (as measured by 
the NO x analyzer), is irradiated, the Og production is affected only 
slightly. Unfortunately, the NO analyzer during these three runs could 
measure oxides of nitrogen as total NO only. It is interesting to note 

A 

that the total NO decreased almost linearly with time. 

A 

Figure 4 shows that when the initial NO is about 65 pphm, the NO 
inhibiting effect on Og production is increased as compared to the re- 
sults with 30 pphm NO (Fig. 3). The Og production delay time is about 
20 minutes. This delay time is defined as being the time duration from 
time zero to a time at which the Og production rate is about 0. 2 pphm 
per minute. The interesting thing here is that this inhibiting effect is 
only temporary, since the total NO in the chamber is again decreasing 
almost linearly with time. 

When the initial NO is 128 pphm (Fig. 5), the NO inhibiting effect 
is increased to the extent that the delay time is now about 80 minutes. 
Again the total NO__ decreases with time, but in a more nonlinear man- 
ner. 


A fourth run, in which the NO was about 60 pphm, was made when 
the NO analyzer was fully operational. These results are shown in 
Fig. 6. The Og reference curve for ”UPC tT air alone was obtained a day 
before this run. These results agree reasonably well with those shown 
in Fig 4, in which the added NO was about 65 pphm. The NO measure- 
ments shown in Fig. 6 can now explain partly why the NO inhibiting ef- 
fect is only temporary. Measurements of NO and NOg indicate that NO 
is rapidly consumed, or converted to NOg (by reaction with Og), and 
that no NO was regenerated. This suggests that for our experimental 
condition, both reaction VI (NOg + h^ — NO + O) and reaction VII (NOg + 

O — NO + Og) were inoperative in recycling NO. The gradual decrease 
in NOg may indicate further reactions involving NOg to form NgO^, di- 
nitrogen pentoxide. This NgOg can, perhaps, be removed by a hetero- 
geneous reaction with trace moisture on the chamber walls to form 
HNOg, nitric acid (Gay and Bufalini (1971)). These reactions will be 
discussed later in more detail. 

Figure 7 is a summary curve showing the effect of varying amounts 
of initial NO on Og production. These experiments indicate that the ob- 
served inhibiting effect of NO is due mostly to gas phase chemical reac- 
tions . 


Effect of Addition of CO, NO, and HgO to 
”UPC” Air on Og Production 

The next series of tests consisted of determining the effect on Og 
production of addition of some initial CO and HgO to "UPC” air contain- 
ing some NO. Prior to these runs involving CO and HgO addition, a run 
was made to determine the effect on Og production of adding initially 
76 pphm NO and about 1200 ppm HgO to air. This result is shown in 
Fig, 8. On this same figure is another reference curve for air alone, 
obtained a day prior to the air + NO + HgO run. As expected, when 
some NO is present, a delay in Og production is observed. Unexpect- 
edly, the 0„ production eventually exceeds that for air alone tempo- 
o 

rarily. The delay time is about 25 minutes, which is about the same 
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as that when the initial NO was 60 or 65 pphm in dry air (Figs. 

4 and 6). 

The effect of initial addition of 50 ppm CO to air containing about 
75 pphm NO and about 1200 ppm HgO is shown in Fig. 9. Two compari- 
son curves are presented, one for ”UPC M air alone and the other for 
M UPC ,f air containing some added NO and HgO. The addition of CO re- 
sults in a slight decrease in Og production delay time to about 21 min- 
utes versus 25 minutes in the absence of CO. No NO measurements 

A 

were obtained for this run due to a breakdown in the NO x analyzer. 

The effect of adding initially 100 ppm CO to TT UPC Tt air containing 
60 pphm NO and about 1200 ppm HgO is shown in Fig. 10. For this run 
the NO x analyzer was fully operational, and NO and NOg measurements 
were made. These results showed that when more CO is added, the 
delay time is now reduced to about 9 minutes versus 21 minutes when 
the added CO was 50 ppm. In addition, the peak Og production was in- 
creased to about 162 pphm. This additional Og production may be at- 
tributed to the action of photosmog type reactions occurring in the reac- 
tion chamber. NO measurements showed that the NO decay is now an 
exponential one, whereas in the absence of CO, the NO decay was more 
linear and slower (see Figs. 3 to 6). Peak NO 2 formation was reached 
in about 25 minutes versus 50 minutes when no CO was present (Fig. 6). 
Again no NO was regenerated just as in the work involving addition of 
NO alone. 

When the CO concentration was increased to about 220 ppm, the ef- 
fect on Og production is shown in Fig. 11. Delay time is now further 
reduced to about 1 minute and Og production of about 172 pphm was 
reached. The NO measurements indicated that both NO and NO de- 

X X 

cays were accelerated, and peak NOg concentration occurred in about 
14 minutes instead of 25 minutes, when 100 ppm CO was added. It was 
observed again that no NO was regenerated. 

Figure 12 is a summary plot showing the effect of additions of vary- 
ing amounts of CO to counteract the NO inhibiting effect on Og produc- 
tion. These experimental results indicate that CO and HgO addition can 
counteract the NO destruction of Og and result in higher Og concentra- 
tion in the chamber. 
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COMPARISON OF EXPERIMENTAL AND 
THEORETICAL RESULTS 
Results With Pure Air 

The steady-state analysis with the simple Chapman mechanism (re- 
actions I-IV) mentioned previously shows that these reactions alone can- 
not account for the observed leveling off of the Og concentration in the 
chamber. An additional destruction reaction for Og has to be used in 
the chemical mechanism. This single reaction is used to represent a 
destruction process which may be a combination of surface and gas phase 
processes. For computations with our general chemical kinetics pro- 
gram, both simple first order and second order kinetics were used for 
the additional destruction reaction to see if either assumption would be 
more consistent with the experimental results. For all computations 
air was assumed to be 79. 05 mole percent Ng and 20. 95 percent Og. 

The exact intensity versus wavelength distribution of our solar 
simulator light source was not known below 250 nm. Thus, none of the 
photochemical rate constants could be computed exactly. Therefore, 
for all computations, reasonable literature values were used at the 
start , When necessary they were changed within reasonable limits to 
give agreement between experimental and computed results. The 
steady-state analysis with the Chapman mechanism shows that reac- 
tion I, the oxygen photolysis, completely controls the rate of formation 
of Og in the early stages of the reaction. Several preliminary computa- 
tions established that any computed Og concentration profile is controlled 
only by two rate constants. These are J^, for the oxygen photolysis, and 
the rate constant for the assumed extra Og destruction reaction in addi- 
tion to reactions III and IV. These last two reactions can in no way ac- 
count for the leveling off of the Og concentration in our experiment. A 
series of computations was done with two different destruction reactions, 
using each separately. The overall destruction reaction used is either 
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or 



(3/2 )0 2 




(d 2 ) 


In performing the computations, the value of was chosen to 

match the initial slope of the Og concentration versus time curve. Then 

the value of or k n was chosen to agree with the final portion of 
U 1 u 2 


the concentration profile. In this analysis a computed curve is chosen 
to match experimental data in a semi quantitative way. It is a curve 
which, in the authors’ judgement, closely represents the trend and mag- 
nitudes of the experimental data. This technique calibrates our system 
and these rate constants can then be used to analyze the experimental 
results in the presence of various amounts of pollutants. Inasmuch as 
there is a small effect of diffusion of ozone from the lighted volume of 
the reaction chamber to the wall, the rate constants found in this pro- 
cedure also include any diffusion effects. 

Computations were also made to match the experimentally observed 
ozone decay in our chamber. This was done to obtain preliminary val- 
ues of k n and k n . Figure 13 shows data for two experiments where 
l>i u 2 


ozone decay was observed. The solar simulator was turned off after 
the ozone concentration reached a steady state value of 120 pphm. The 
mechanism used to compute the expected decay curve consists of reac- 
tions II, IV, and either or D 2 . The value of k^ or k^ was found 

to completely control the ozone decay curve. The two curves in Fig. 13 
show the best results obtained using either reaction or D 2 . The 

values obtained are k. 


n = 4x10 ^ sec" 1 and k^ = 5x10^ cm 3 /mole- 
l>i d 2 


sec. The curve using reaction (uni molecular decomposition) follows 

the trend of the experimental data somewhat better than the curve using 
the bi molecular destruction, reaction D 2 « 
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The ozone destruction mechanism may not be exactly the same with 
the light on as with it off. Therefore, both reactions and Dg were 
used (separately) in two series of computations to match the ozone for- 
mation experimental results. The value of kp. or kp. was first set 

U 1 u 2 


at the value given above. Then the oxygen photolysis rate constant 
was changed to obtain good agreement between experimental and com- 
puted results for the beginning of the reaction. The k n or kp. value 

U 1 u 2 


was then varied slightly to match the experimental data for the end of 
the reaction. Figure 14 shows the results of these trials. The data 
points are for three different ozone formation experiments performed 
over a 2 -month period before the chamber surface was exposed to any 
pollutant. The two curves are the computed curves obtained using the 
two different destruction reactions, and Dg. Rate constants used 


-9 -1 -4 -1 

for these computations are = 1x10 sec , k^ = 3x10 sec , 

c o 1 

and k n = 4x10 cm /mole- sec. The two curves are very much the 
u 2 


same. Either reproduces the experimental data quite well. These 
computations do not distinguish whether a first or second order Og 
destruction reaction occurs in the reaction chamber. Therefore, for 
the computations in the presence of pollutants, both reactions were con- 
sidered. All computations were performed first using reaction and 
then repeated using reaction Dg. The rate constants used for all reac- 
tions considered in this work are listed in Table I. 


Results With Nitric Oxide 

Computations were performed for the case of nitric oxide added to 
dry air. The first mechanism used consists of reactions I through XI 
plus reaction Dg. The computed results did not agree qualitatively or 
quantitatively with the experimental results. The ozone concentration 
leveled off at a value much lower than that observed experimentally. No 
rate variations within estimated uncertainties changed these qualitative 
results. Adding two additional reactions, however, had a significant 
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effect on the computed ozone profile. These reactions, which involve 
nitrogen pentoxide, NgOg, are 

^22 

N0 2 + NOg + M > N 2 0 5 + M (XXH) 


^ 23 a 

hv + N 2 0 5 — > 2N0 2 + O (XXIIIa) 


J 23b 

> no 2 + NOg (xxnib) 

A comparison of computed and experimental results using reactions 
XXII and XXIIIa is shown in Fig. 15(a) and (b). Part (a) shows computa- 
tions using destruction reaction Dj and part (b) shows computations using 
reaction D 2 . For either reaction, the computed curves follow the quali- 
tative trends of the experimental curves. For 30 pphm of NO, the com- 
puted delay in 0 3 formation is longer than the observed delay. For the 
higher initial NO concentrations, the computed delay is shorter than 
that observed. 

Analysis of the rates of formation and destruction of ozone by indi- 
vidual reactions shows that different ones are important at different 
times. At the beginning of the reaction the 0 3 formed by reaction II is 
destroyed primarily by the reaction 

NO + O s - N0 2 + 0 2 (V) 

As ozone is built up it is also destroyed by 

h^+ 0 3 - 0 + 0 2 (III) 

Another destruction reaction that is relatively unimportant compared to 
the previous two reactions is 

N0 2 + 0 3 NOg + 0 2 


(VIII) 
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The additional destruction reactions and Dg are even less important 
than reaction VIII in destroying Og at the start of the reaction. 

At the end of the reaction, Og concentration levels off because reac- 
tion II is balanced primarily by the effect of the photolysis destruction, 
reaction III. The destruction reaction or Dg has a much smaller ef- 
fect in determining the final Og level for this system than for the pure 
air system. These computations support the previous interpretation of 
the experimental results. It was stated previously that the experimental 
results should be representative of homogeneous gas phase reaction in 
the chamber. The effect of the chamber walls is small during most of 
the reaction time. 

Figure 16 shows several computed species concentration histories 
for the air + 65 pphm nitric oxide case. Comparison with the experi- 
mental curves shown in Fig. 6 indicates good qualitative agreement in 
general. Nitric oxide is rapidly converted to NOg. Nitrogen dioxide 
concentration peaks at about the maximum slope of the ozone concentra- 
tion curve, and then decreases. The experimental decrease of NC >2 
(Fig. 6) is slower than that computed because the NOg measurement 
probably includes higher oxides of nitrogen. Computations show that 
NgO^ is formed slowly after about 20 minutes and increases to a steady- 
state value. 

The preceding discussion refers to computations using reaction 
XXIIIa, the NgO^ photolysis reaction giving the products 2NOg + O. All 
results are essentially the same when the computations are repeated 
using reaction XXIIIb, which gives the products NOg + NOg. The com- 
puted concentration profiles were compared with those of the early com- 
putations, without any NgO^ reactions present. It was observed that 
NOg concentration remains high during most of the reaction when it can- 
not react to form NgO^. If NOg concentration stays high, computations 
show that it will destroy Og by the direct NOg + Og reaction. The present 
analysis gives no information about the products of the NgO^ photolysis. 
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Results With NO, CO, and H 2 0 

The mechanism for Og chemistry in the presence of CO, NO, and 
HgO includes reactions I to XXIII plus reactions which provide OH radi- 
cals, The latter are necessary for starting the chain (reactions XII 
to XV) which converts NO to N0 2 without the destruction of Og. The 
following photolysis reactions were used as the source for OH: 



d 24 



hv + H 2 0 

— ^ 

OH + H 

(XXIV) 


J 25 



hv + Hg Og 


OH + OH 

(XXV) 


J 26 



hi/ + HN0 2 

> 

NO + OH 

(XXVI) 


J 27 



h v + HOp 

> 

0 + OH 

(XXVII) 


Hydrogen peroxide, HpOg, is formed by the reaction 

k 28 

ho 2 + ho 2 — > h 2 o 2 + 0 2 (XXVIII) 

Rate constants for the photolysis reactions were chosen to give 
satisfactory agreement with the experimental results. The values used 
are generally within a factor of 4 of those given by McConnell and 
McElroy (1973) for an altitude of 80 km. However, J^, the HgO 
photolysis rate constant, was taken as l/10th of their value. Computa- 
tions show that the two important sources of OH in our mechanism are 
the H 2 0 and H 2 0 2 photolyses. 

The rate constants used in the dry air plus NO computations were 
not changed for the computations with CO and HgO also present. Com- 
puted results were found to be only slightly different using Og destruction 
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reaction from results using reaction Dg. Inasmuch as the 10 per- 
cent higher final Og concentration using Dg was more in agreement with 
experiment, the computations to be discussed use reaction Dg. 

Figure 17 shows both the experimental and typical computed results 
for Og production in air containing 50 ppm CO, 75 pphm NO, and 1200 
ppm ^O. Also shown for comparison are experimental and computed 
Og concentration curves for TT UPC’ T air + 65 pphm NO added initially. 
These last two curves are the same ones shown in Fig. 15(b). There is 
poor quantitative agreement between experimental and computed curves 
in the presence of CO and I^O. However, the computed results show 
the same qualitative effect of CO that is observed experimentally. The 
computed delay in Og formation is much less in the presence of CO and 
H 2 O. It is clear that the computations do not fully explain the chemical 
reaction in the presence of CO and HgO. The slope of the main portion 
of the computed curve is much lower than that for the experimental 
curve. Moreover, the steady state Og concentration is also lower on the 
computed curve. Another problem is the failure of the computations to 
show increased effectiveness of higher CO initial concentrations. The 
computed results for 200 ppm CO present show a significantly smaller 
difference from the results for 50 ppm CO than is observed experimen- 
tally. The qualitative trend is shown, but there is no quantitative agree- 
ment. Many rate constant variations showed that the actual delay in O^ 
formation is quite sensitive to the rate constants and for photo- 
chemical formation of OH radical. However, for all perturbations of 
other rate constants within their estimated uncertainty, the slope of the 
main portion of the O^ concentration curve and the final O^ level were 
changed only slightly. 

One final point concerning the reactions of NO 2 formed from the 
original NO, should be mentioned. We have discussed earlier the pos- 
sibility of the removal of higher oxides of nitrogen from the gas phase 
as HNOg condensed on the chamber walls (Gay and Bufalini (1971)). 

This idea provides one possible reason for the disagreement between 
computed and experimental Og production. The computations do not 
allow the removal of higher oxides of nitrogen from the gas phase. 
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Our reaction mechanism gives only a relatively small amount of indirect 
conversion of NgO,- to HNOg via reactions XXIII and XX. This is cer- 
tainly unrealistic. Our gas phase mechanism does not properly explain 
the subsequent reactions of NOg. Additional computations must be done 
to better explain our experimental results. 


CONCLUDING REMARKS 

The main concern about the effect of nitric oxide emission in the 
stratosphere has been its possible long lasting destruction of ozone. 

Many papers (e.g., Johnston (1971), Westenberg (1972)) have predicted 
that NO destroys ozone catalytieally, i. e. , without itself being destroyed. 
These conclusions were based on approximate steady state analyses and 
the assumption that all the NOg formed by reaction V would be photolyzed 
back to NO by reaction VI. The present laboratory scale experiment in- 
dicates that NO is destroyed while destroying Og. The Og formation 
reaches its usual level (in the absence of pollutants) soon after most of 
the NO has been converted to NOg and other oxides of nitrogen. Even 
though NOg is partially photolyzed back to NO, this reaction also pro- 
duces oxygen atoms. These, along with the O atoms from the photolysis 
of Og, form more ozone than the regenerated NO can destroy. The com- 
putations for dry air with NO added show that NO is destroyed much 
faster by the NO + Og reaction than it is regenerated by NOg photolysis 
in the early stages of the reaction. In the latter part of the reaction NO 
concentration becomes quite low. Then the NO concentration stays at 
this low level because its formation by NOg photolysis (and the reaction 
O + NOg — NO + Og) is exactly balanced by its destruction in the NO + Og 
reaction. At this time the Og forming reactions cause the Og level to 
build up to its usual level in the absence of NO. 

The effect of CO and HgO, which increases Og concentration in the 
present experiment, is not completely understood. Since these products 
are present, along with NO, in engine exhausts, the matter of equilib- 
rium Og concentration in the presence of engine discharges is still un- 
certain. 
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An important question is whether our laboratory results can be ap- 
plied to low pressure conditions. To partly answer this question, the 
experiments and the theoretical computations just described will be re- 
peated at reduced pressures, down to about 0.3 atm. The latter pres- 
sure is the expected lower limit, based on preliminary tests of our 
modified Og analyzer. If the theoretical modelling is successful at this 
reduced pressure, the same chemical mechanism will be used for theo- 
retical computations at stratospheric pressures. The attempt will be 
made to predict the trend of Og formation in the presence of pollutants 
down to 0.05 atm. The ratios of the various pollutants will be adjusted 
to simulate upper atmospheric concentrations. 
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Figure 1. - Simulator and reaction chamber. 



Figure 3. - Ozone production for "UPC" air containing 
30 pphm NO. 



Figure 2. - Ozone production for "UPC" air. 
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Figure 5. - Ozone production for "UPC" air containing 128 pphm NO. 
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Figure 6 . - 0 3 production for "UPC" air containing 
60 pphm NO. 
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Figure 7. - O 3 production for "UPC" air containing 
various initial amounts of NO. 
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Figure 8. - Ozone production for "UPC" air containing 
76 pphm NO and 1200 ppm H 2 0. 
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Figure 9. - Ozone production for "UPC" air containing about 
50 ppm CO, 75 pphm NO, and 1200 ppm H 2 O. 



Figure 10. - Ozone production for "UPC" air containing about 
100 ppm CO, 60 pphm NO, and 1200 ppm HjA 



Figure 11. - Ozone production for "UPC" air containing 
about 220 ppm CO, 65 pphm NO, and 1200 ppm F^O. 
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Figure 12. - Ozone production when various amounts of CO 
are added to "UPC" air containing some NO and H 2 0. 
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Figure 13. - Ozone destruction in "UPC" air; experimental and 
theoretical. 
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Figure 14. - Ozone formation in "UPC" air ; experimental and 
theoretical. 
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Figure 15. - Comparison of experimental and theoretical 
ozone formation in air + nitric oxide mixtures. 
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Figure 16. - Species concentration vs. time for irradiation of 
air + 65 pphm nitric oxide; computed results using uni- 
molecular destruction reaction, O 3 - (3/2) 0 2 - 
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TABLE I. - CHEMICAL REACTIONS AND RATE CONSTANT! 


Number 

Reaction 

Preexponential 
factor, * 

A 


hv + 0 -? -* 0 + 0 
0 + 0 2 + 1 M “• 0? + M 
hv + — 0 + Og 

0 + O 3 ■* o 2 1 - o 2 

1 x 10 ' 9 
3.83X10 13 ■ 
3.52x10 ~ 6 
l.lxlO 13 

V 

NO + 0 3 ** N0 2 + 0 2 

5. 5x10 11 

VI 

hv + N0-? - NO + 0 

1 x 10 ” 3 

VII 

0 + NQ 2 NO + 0 2 

1 ? 

5.5x10 

VIII 

i\ 0 o + o 3 rro 3 + 

3x10 7 

IX 

2N0 + 0q - 2N0 2 

1.2X10 9 

X 

NO + 0 + M - N0 2 + M 

9.4xl0 14 

XI 

hv + N0 3 - NO + 0 2 

1 x 10 ' 3 

XII 

CO -r OH - C0 2 + H 

8 . 7x10 10 

XIII 

H t 0-, + M -* H0 o + M 

1.59x10 15 

XIV 

NO + H0 2 -* N0-? + OH 

3.6X10 11 

XV 

H -r N0 2 - NO -r OH 

2 . 9X10 13 

XVI 

0*1 1- 0 3 -* HO 9 ■+* 0 ^ 

3X10 11 

XVII 

H0-? t Oj ■* OH + 20 2 

3x10 10 

XVIII 

NO + OH +■ H -» HNOg + M 

3xl0 16 

XIX 

C*i ■*“ *lu0o H^Q 

8.4X10 11 

XX 

N0g OH + M *♦ HNO 3 + M 

8 xl 0 16 

XXI 

OH + M;C 3 - HoO * NO 3 

8.4X10 11 

XXII 

N0 2 + N0 3 i- M - N 2 0 s + M 

1 x 10 18 

XXI I la 

hv + 2 N 0 2 + 0 

1 x 10 ’ 4 

XXI I lb 

hv t - ..p 0 g ”* N0-? + N0g 

7. 7x10 " 5 

mv 

hv + n 2 0 -* H + OH 

1 x 10 " 7 

XXV 

hv + H o 0 o -* OH + OH 

5. 7x10 " 5 

XXVI 

hv + a.'0 2 - NO + OH 

1.3xl0 ‘ 3 

XXVII 

hv + H0 2 -* 0 + OH 

3. 9x10 " 4 

XXVIII 

MU 2 "f iiOp HgOg ■** 

5xl0 12 


exponent, 

n 


energy, 

kcal/mole 



See text 

Garvin ana Hampton (1074) M = 
Westenberg (1972)7 

Garvin and Hampton (1074) 
Garvin and Hampton (1074) 
McConnell and HcElroy (1973) + 
Garvin and Hampton (1074;** 
Garvin and Hampton (1974)** 
Garvin and Hampton (1074) 
Schofield (1067) 

McConnell and HcElroy (1573 )"*" 
Garvin and Hampton (1074)** 
Baulch, et al. (i960) H = Ar* 
Davis (1972)7,-** 

Gam/in and Hampton (1974) ** 
Garvin and Hampton (I574) t 
Garvin and Hampton (1974)7 
Garvin and Hampton (1074)'"'* 
Hampson (1972) 

Simonaitis and Heicklen (1972)' 
Hampson (1972) 

Garvin (1973)** 

Johnston (lOVl ) 7 
McConnell and HcElroy (1973) 
McConnell and HcElroy (1973)* 
McConnell and HcElroy (1973) 
McConnell and HcElroy (1973) 
McConnell and HcElroy (1973) 
Garvin and Hampson (1974)7 


"^Constants in the equation k = AT a exp(.-E a /RT) * Units of k are sec" x for photochemical reaction, 
cm 3 /mole sec for binolecular reaction; and cm°/mole 2 sec for termolecular reaction. 

* Adjusted from value given in reference, 

*Third body relative efficiencies, N 2 /Ar and 0 2 /Ar, for this reaction are taken as 2. A more recent 
recommended equation (Garvin and Hampson (1974)) gives a value about 20 percent lower at T = 300 H. 
"^Temperature = 298 K. 









